home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / fractal / kaos.lha / modellib / pp2harvest_def.c < prev    next >
Encoding:
Text File  |  1990-01-16  |  1.1 KB  |  59 lines
  1.  
  2. int pp2harvest_init()
  3. {
  4.     title_label = "2-Species Predator Prey with Harvesting";
  5.  
  6.     mapping_on = 0;
  7.     inverse_on = 1;
  8.     fderiv_on = 0;
  9.     enable_polar = 0;
  10.     enable_period = 0;
  11.  
  12.     var_dim = 2;
  13.     param_dim = 4;
  14.     func_dim = 2;
  15.  
  16.     (void) malloc_init();
  17.  
  18.     var_label[0] = "u1";
  19.     var_label[1] = "u2";
  20.     param_label[0] = "p1";
  21.     param_label[1] = "p2";
  22.     param_label[2] = "p3";
  23.     param_label[3] = "p4";
  24.     func_label[0] = "t";
  25.     func_label[1] = "Undefined";
  26.  
  27.     param[0] = 2;
  28.     param[1] = 1;
  29.     param[2] = 0;
  30.     param[3] = 0;
  31.     var_i[0] = 0;
  32.     var_i[1] = 0;
  33.  
  34.     param_min[0]= -5; param_max[0]= 5;
  35.     param_min[1]= -5; param_max[1]= 5;
  36.     param_min[2]= -5; param_max[2]= 5;
  37.     param_min[3]= -5; param_max[3]= 5;
  38.     var_min[0]= -5; var_max[0]= 5;
  39.     var_min[1]= -5; var_max[1]= 5;
  40.     func_min[0]= -5; func_max[0]= 5;
  41.  
  42.     f_p = pp2harvest_f;
  43.     func_p = pp2harvest_func;
  44. }
  45. int pp2harvest_f(f,index,x,p,t,dim)
  46. int index,dim;
  47. double f[],x[],p[],t;
  48. {
  49.     double exp();
  50.     f[0] = p[1]*x[0]*(1-x[0]) - x[0]*x[1] - p[0]*(1-exp(-p[2]*x[0]));
  51.     f[1] = -x[1] + p[3]*x[0]*x[1];
  52. }
  53. int pp2harvest_func(f,x,p,t,dim)
  54. double f[],x[],p[],t;
  55. int dim;
  56. {
  57.     f[0] = t;
  58. }
  59.